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ABSTRACT: 


The report presents a computer program for prediction of performance 
of single-stage axial turbines of given geometry. The three-dimensional 
method developed by Vavra is applied, taking account of streamline 
curvatures and slopes, as well as enthalpy and entropy gradients in 
the solutions of the equation of motion, and of boundary layer thick- 
nesses in the continuity equation. 


A choice among five different loss correlation methods and two 
flow angles correlations is offered. Loss coefficients and flow 
angles are automatically calculated from blading geometry and actual 
flow conditions for every streamline, according to the selected 
correlation method. 


A fair agreement of predictions with several actual turbines 
experimental results was found in ref. Ley where also the applicability 
of different available correlations is discussed. 
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LIST OF SYMBOLS 


Area (in?) 

Tip clearance area (in@) 

Throat opening of blade channel (in) 
Non-dimensional throat opening, a/ a, 

Blade chord, see Fig. 6 (in) 

Rotor tip clearance (in) 

Specific heat at constant pressure (BTU/1b,-°R) 
Infinitesimal area, perpendicular to the flow (in?) 
Blade radius of curvature, see Fig. 6 (in) 
Frictional force per unit mass (1b/ slug) 

Total enthalpy (BTU/1bm) 

Boundary layer form factor 

Boundary layer energy factor 

Coefficient, see eq. (31) 

Curvature coefficient, see eq. (16) 

Unit vector 

Incidence angle 

Mean distance between the stations O and 1 and 1 and 2 (in) 
Mach number 

Boundary layer exponent 

Mass flow rate (slug/sec) 

Mass flow rate passing thru the rotor tip clearance (slug/sec) 


Number of conditions for which performance values have to be 
found 


Blade opening (in) 

Total pressure (psi) 

Statice pressure (psi) 

Heat exchanged with external sources per unit mass (BTU/1b) 
Gas constant (ft-1b/1b,-°R) 

Radius at blade hub (in) 

Radius at blade tip (in) 

Radius (in) 

Entropy (BTU/1b-°R) 


Non-dimensional entropy, S/¢p 


= 
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Blade maximum thickness, see Fig. 6 (in) 
Static temperature (°R) 
Total temperature (°R) 


Thickness of blade trailing edge, in peripheral direction, 
see Fig. 6 (in) 


Thickness of blade trailing edge, in normal direction, 
see Fig. 6 (in) 


Time (sec. ) 

Peripheral velocity (ft/sec) 

Absolute velocity (ft/sec) 

Relative velocity (ft/sec) 

Non-dimensional radius, r/rp 

See eq. (30) 

Non-dimensional axial velocity, Va/Van 

Total pressure loss coefficient, see eq. (40 ) 
Number of blades 
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GREEK SYMBOLS 


Absolute gas flow angle, see Fig. 5 (deg or radians) 
Relative gas flow angle, see Fig. 5 (deg or radians) 
Pressure ratio, see eq. (44) 

Blade inlet angle, see Fig. 5 (deg or radians) 


Specific heat ratio 


Decrease in efficiency, due to rotor tip leakage, see eq. (47) 


Streamline displacement, see Fig. 1 
Streamline displacement, see Fig. 1 
Kinetic energy loss coefficient 


Kinetic energy loss coefficient used in the continuity 
equation, see eq. (28) 


Kinetic energy loss coefficient for low Mach numbers, 
see eq. (42) 


Kinetic energy loss coefficient due to rotor tip leakage, 


see eq. (47) 

Total/static efficiency 

Total/total efficiency 

Angle between two radial planes (deg. or radians) 


Angle between the flow and the axis of the turbine in a 
meridional plane (deg. or radians) 


Gas viscosity (1b/sec-ft) 

Area restriction factor, see eq. (27) 
Gas density (1b,,/ ft?) 
Non-dimensional flow function 


Coefficient in Traupel's method of loss correlation, 


see eq. (46) 
Angular velocity (radians/sec) 
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SUBSCRIPTS 


Axial 

Equivalent 

Isentropic expansion from total inlet conditions 
Mean streamline 

Meridional 

Rotor 

Radial direction, cylindrical coordinates 
Stator 

Axial direction, cylindrical coordinates 
Peripheral direction, cylindrical coordinates 
Station ahead the stator 

Station between stator and rotor 


Station after the rotor 


INTRODUCTION 


An accurate knowledge of performance achievable with turbines 
of fixed geometry, both at design and off-design conditions, is 
important in advanced propulsion and power systems. The more accurate 
the prevision is, the more meaningful the optimization of the turbine 
itself and of the overall system will be. 

A first way to solve this problem consists in the testing of a 
prototype, but it presents several disadvantages, namely it is 
expensive and time consuming. Furthermore, if effects of turbine 
geometry changes have to be investigated, different prototypes 
will be required. 

An alternative solution to the problem consists in theoretical 
methods of prediction. The present availability of high speed 
computers allows one to solve systems of differential equations and 
to perform a large number of iterations in a short time. Moreover 
computer programs give the possibility of changing imput data,so that 
a large variety of different turbine configurations and of working 
conditions can be investigated. 

This report presents a computer program, derived by the 
three-dimensional calculation method given first by Vavra [1] and 
then developed by Eckert [2] and Harrison Lewes This method takes 
account of enthalpy and entropy gradients and of streamline curvatures 
and slopes in equations of motion, and of boundary layer thicknesses 
in continuity equation. It is applicable to single-stage axial 
turbines with unchocked flow. 

A critical point in all theoretical performance calculations 
is the prediction of the losses that occur in the blade rows. In 
the present program, one can choose among five different loss correla- 
tion methods. For each method pertinent subroutines compute loss 
coefficients from blading geometry and flow characteristics. A 
comparative study of these methods can be found in [4], where they 
were applied to actual turbines of different design practices, in a 


wide range of all significant flow and geometrical properties. 


The comparisons presented in [4] with experimental performance 


showed a fair agreement of predictions with all turbine stages 
that were examined. 

The present report is divided into four sections: in the first 
the used equations are set up and necessary assumptions are discussed; 
in the second the method of solution is described; in the third 
different losses and discharge flow angles available prediction methods 
are described; and in the fourth the instructions for the use of 
computer program are given, with an example of imput and output data. 
Eventually in the appendix the FORTRAN list of the MAIN routine and 
of the 28 subroutines is given. 

METHOD OF ANALYSIS 

The flow in turbomachines is ruled by the three fundamental 
laws of fluid dynamics, namely: 

Equation of Motion 

Energy Equation 

Equation of Continuity. 


All proper three-dimensional calculation methods proposed 
in the technical literature solve the differential equations derived 


by these laws. Differences among these methods are due to the three 
following reasons: 
- different stations are chosen where equations are solved; 
- different hypotheses are formulated to describe the real 
processes occurring in the machine; 
- different numerical methods are applied in the solutions of 
the equations. 
As far as the first problem is concerned in the present method, 


the equation of motion is applied at three stations, namely ahead of 


the stator, after the stator, and after the rotor, indicated respectively 


with subscripts 0, 1, 2 in Fig. 1. On the contrary, since blade geometry 


is given, the continuity equation is applied at blade throat sections 


1* and 2* in Fig. 1, as suggested in [5]. 


146) 


those occurring in the boundary layers along the blade surfaces, and 
flow frictional losses between blades can therefore be neglected. 

Assumption 4) is reasonable for uncooled blades, where energy 
changes due to heat transfer are negligible with respect to those 
related to velocity changes. If blades are cooled, the additional 
term q in the energy equation must be specified. Assumption 5) is 
valid in most cases. If the thermodynamic processes occur in ranges 
where real gas effects are not too strong, average values of gas 
properties give good approximations. 

Some additional assumptions, regarding streamlines slope and 
curvature, are discussed in the section pertaining to the equation 
of motion. 

As far as the third problem is concerned, the numerical 
small difference method was found to give a good convergence, and 
was therefore applied for the solution of the differential equations. 

The following sections give the derivations of the equations 
with the above-listed assumptions. 


Equation of motion 


The general equation of motion for relative flow in vectorial forn, 


is: 

Set Hp = WX(VxW+ 20)+TVS+ fp (1) 
TPO (2) 

For absolute flow, where ve=W (3)  , equation (1) 
H =H, (4) 

becomes: 

QU +VH=Vx(VxV)+T S84 fe (1a) 

Wi r, hoe Me CN i 

ith hypothesis 1) SOs St TO 
and hypothesis 2) ae = 0 


ae: 


This allows one to more accurately evaluate the pressure drops 
occurring in blade channels, since in this way it is possible to 
account for boundary layer thicknesses. This problem is particularly 
important at high Mach numbers where a small change in flow area causes 
large differences in pressure ratios. 

As far as the second problem is concerned, the following 
assumptions are made: 

1) Steady flow 

2) Frictionless flow at the stations where the equations of 
motion are solved. 

3) Axisymmetric flow at the stations where the equations of 
motion are solved. 

4) Adiabatic flow. 

5) The working fluid is a perfect gas, namely a gas having 


the equation of state: 


p=PRT 


and specific heat at constant pressure is independent of temperature 
and pressure. 

Hypotheses 1) and 3) are interrelated since the flow can be steady 
in both rotor and stator only if it is axisymmetric. Although they 
are correct only for cascades with an infinite number of blades, 
axisymmetry and steady conditions are necessary assumptions since any 
other hypothesis would require unknown quantitative information on 
flow inside the blade rows. One problem connected with these assumptions 
is that downstream conditions cannot affect the flow upstream. There- 
fore pressure distributions at rotor exit can deviate from those imposed 
by the discharge conditions. 

With assumption 2), all entropy changes are supposed to occur in 
the blade rows ahead of the considered station. This assumption can 
be justified by the fact that velocity gradients in the regions, where 


the equation of motion is applied, are in all probability smaller than 


¥2 


equations (1) and (la) become respectively: 


VHp=Wx(VxW+ 20)+TVS (5) 


= 


VH=Vx(VxV)+TVS (5a) 


The following derivations for equation (5) are valid also for 
equation (5a), if substitutions (2), (3), (4) are performed. 

It is convenient to express equation (5) in cylindrical coordinates 
(see Fig. 1 for symbols). Equating the three components of equation (5), 


the following equations are obtained: 


a. tangential component: 


ee a ee ee | ay TES 
mao. rr = 08 “oz Bley, air ~ 68 a Gr | £00 
b. meridional component: 
a fl (7) 
a) 6 he Om or A Dow me. Oz Pe ay, 
ec. radial component: 
OH, Wy | a(rWy) dW | OW, Wyss 3 
rote che SO me pe eee 8S 
or  r L ae Bo L Oz 9 oR | ee Sr (8) 
with hypothesis 3) e (Co) =O 
yp RG) = 

ene bypothesis~ 4) Syd -t°=°V Hp 240 
the system of equations (6), (7), (8) reduces to equation (9): 
oH Ww. 9(rW_) OW. QW 3 

R UL UL r A Js 
mec re = VA Se eo ot aed Eee (9) 
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Using substitution: Hp = Hp + aa (10) 


with hypotheses 5) H = Cp T and 4) H, = const (see Fig. 2) 


== = Zou (em) 


and introducing the non-dimensional quantities: 


W 
Yoo (12) 
Am 
r 
= (13) 
Tn 
ox = = (14) 
Pp 
Equation (9) becomes: 
aye er, oW. 
CRG 2 m or 1 ag* op 
re + cos 6 (- Tp es a sa) + 2 tan Pp =, 
z 
ae 4 sin 6B cos B aU oe Pp 2 aoe p Hp 
ar RE Fei, eaee Lhtebigetge winne Gee eiee: aameaal 
Am Wa Ys Wam 
—-2 H cos‘pB Ke 
+ — - sin6 as 0 (15) 
Yet ye ox 
Am 


OW: 
In this equation, terms containing cos 4’ and > depend respectively 


on streamlines slope and curvature. Vavra [1] showed that the most 


reasonable assumptions, when only three points of a streamline are known, 


are the following (see Fig. 1 for symbols): 


SW 6 
os Wie Ce ia (16 ) 
Wa Z it 
2 
2 = L 
NS ae ee (17) 
Le cs =) 
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In equation (16) the positive sign holds for station (2) and the negative 
one for station (1), for 6r positive as shown in Fig. 1. Suggested values 
for k are between 4 and 6. Finally with equations (16) arid C7). ntor (15 ) 


the following equation is obtained: 


d (1n Y°) r tee Gy ase 
cee | _ cose8 | 2k er =. | 
ax m Tf? i ax 
dp 2 4 Un sin8 cosf 
- 2 tan 6 = - sin p - rien = 
m0! U cos“B -e B cos“6 
_ ™m 2 cos*B dHp i mbar) asx | 
ee roe heer) 
Wy Y Wan, Yo ax i eee oe 
Am Am Am 


For the station after the stator row, the following equation can be 


derived with (2), (3) and (4) in equation (18): 


cos 
7 Li Le dx 
2 tan & oe eee og | 200s" @ jel 2 H cos< a 
c- 1 pais bs — a 2 dx ~ 5 - 
ye Vien ye gee 
% 
_ a. = (18a) 


Equation(18a) could also have been derived directly from equation (5a). 
In equations (18) and (18a) the last term is depending on flow 
irreversibilities due to frictional losses. The specific entropy S* will 


be calculated, after Lit by: 


2 ee 
Y Wan 


2] 
Pai COSs 
Oe Sra (19) 
Y Wam 
1- ———__-—— 
2 Hp cos B (1-<) 
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=. Vi 
sS* = ln | ————— + § (19a) 


2 Ho cos<8(1-£) 


where € represents the kinetic energy loss coefficient which is 
established by the chosen loss correlation method. 
Energy equation 

The energy equation in a relative flow is given by: 


ie 


Wat-V (Hp) = Wdt-Y(q) - 3¢ (5) at (20) 
and in an absolute flow, with (2), (3) and (4) in (20), : 

3 ve 
Vat-V (Hp) = Vdt-V(q) - 4 (3) dt (20a) 


where q is the heat added to a unit mass particle from external 
sources. With assumptions 1) a/ct( ) = 0 and 4) q = 0, equations 
(20) and (20a) become: 


Wat-V (Hp) = 0 (2) 


Vdar-V (H) = O (21a) 


Equations (21) and (2la) indicate that for steady, adiabatic flow 
the relative total enthalpy is constant along any streamline. This 
does not mean that V (Hp) = O, since relative total enthalpy can assume 
different values for different streamlines. 
Equation of continuity 

Pressure ratios across a stage necessary to discharge a given 
flow rate are function of inlet Mach number, blade openings and the 
losses occurring in the blades. For an isentropic flow the mass flow 


rate passing thru a section dAn perpendicular to the flow is given by: 
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Om, = Pig Mpprahs (22) 


where Pig and ve are respectively the density and the velocity 


S 
corresponding to an isentropic process. It can be shown that 


equation (22) can be expressed for a perfect gas (see assumption 5) 


by: 

qa, = 2 aA, (23) 
pie 

where Pp, Tp are, respectively, total pressure and total temperature 


(constant for isoentropic and adiabatic flows) and 4 is a function of 


pressure ratio and gas isentropic exponent, namely: 


(2h) 


With the symbols in Fig. 3, assuming the pressure to be constant along 


the throat "a" equation (23) becomes: 


diis = ——— facosrAdr (25) 


Integrating equation (25) along the radius and for z blades the 
flow rate passing thru the stage is: 


é a cos \ dr (26 ) 


In actuality, the flow rate will be less than m;,, because of the losses 
occurring in the blade channels. 
A so-called restriction factor: 
d wm 
SS 


(27) 


is now introduced. 


Ay 


In [5] the following relationship between & and the kinetic 


energy loss coefficient (* is derived: 


g 21 ee (28) 


where H and K are respectively the boundary layer wake form factor 
and the energy factor. For boundary layer profiles following the 


5 F ; 
power law, the ratio; can be expressed in power series: 


[eo) 
= 


12 (lex) vb xe Tas(2i-1)m| 


si E (29) 
(oo) i aa eo alae ‘ a 
(1- Xe ) = | Xe" aia (1-Xe) oi | xe 7 /1+(2i+1)m | 
where: 
Y-1/ y 
Xe = 1 - (P/Pp) (30) 


and m is the exponent of the power law for the boundary layer profile 
(for incompressible flow,m depends on Reynolds number and varies between 
1/7 and 1/10). Experience has shown that the loss coefficient ¢* has 

to be taken less than the value Srop that corresponds to the overall 
energy losses that are given by the correlation methods. In the present 
program one of the three following different assumptions for ¢* can be 


chosen: 1) C* = Snonp 5 2) C¥ = SCmop 3 3) &* = CpROFTIE - 


While assumption 1) is certainly conservative, it is difficult to state 
which of assumptions 2) or 3) is the best one. Good agreement between 
experimental values and predicted flow rates for all considered turbines 
at low tip clearance values was found in ny where assumption 2) was 
adopted. Assumption 3) was however found in closer agreement with 


test data for large values of rotor tip clearances. 
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Tip leakage flow rate will be calculated by using the relation: 


dt Po 
Mar = (< 6 cos \) A (31. 
SL io. = TIP “CL ) 
where: Ajgy = 27 R, cx (32) 


Since the flow is not perpendicular to the leakage area, and since the 
restriction factor Srp does not account for boundary layer blockage 
occurring on shroud, the leakage flow rate must be divided by a 
coefficient Kor, (larger than 1), assumed equal to 2 in the present 
program. 

With equations (29) and (33), the overall flow rate passing thru 


a blade row is given by: 


i : (3) 
- = : ' . ak 1 ci 2 
lipop = Mh + the, = Z 5 Pp? Peon ae Se cosh) 7p lee 
V 
Ry oe CL R Tp 


Xn (34) 
h = z Ry a, : Ce ia SEO ¢ @cos \). dx + ews nicl ( “ $ cosh) 
- ————— x "TIP x Ss Pe 
Xy Vv Rg Tn Koy, xX, 2a ‘Rs Ee 
where a = = (35) 
m 


METHOD OF SOLUTION 


Since no direct solution of governing equations is possible, an 
iterative method is followed. At first, solutions are found without 
considering the influence of streamline curvature and slope. The 
derived radial shifts are then used for the second and final cycles of 
calculation. In Fig. 4 the scheme followed in the computer program is 
indicated. The names in brekets correspond to the names of the 


subroutines which perform the operations described next to the brackets. 
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The flow ahead of the stator is supposed to be uniform, so that the radii of 
streamlines dividing the flow rate in equal parts can be easily computed. 


For stator and rotor rows the following functions are supposed to be 


known: 

STATOR ROTOR 

a = Oy (r, flow condition) Bo = Bo (r, flow condition) 
a, = a) (r) A = Ay (r) 

Ga os (r, flow condition) CR = Sp (rv, flow condition) 


All derivatives used in the following analysis are computed by 
finite difference methods. The computer program calculates flow 
conditions along 5 streamlines. However the same scheme of calcu- 
lation can be applied to any number of streamlines. By assuming the 
inlet Mach number M,, the flow rate passing thru section O can easily be 
determined. Radial position of streamlines at station 1, R,, and the axial 
velocity at mean streamline, Vyjj,,. are also assumed a priori, and will 
be checked later. 
For any given values of V , AR, and 6r, equation (18a) can be 


solved by iteration. In fact it can be reduced to: 


a(n Y°) 
Ee = dale.) (36) 
F355 =) 7 

Equation (36) integrated, with boundary condition|y = 1 (37) 
gives: . = 

in Y= | 4 ax (38) 

al 
and 
pee 
Y = qe ato (39) 
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Solutions are found by trial and error, assuming approximate values 
of Y(x) and verifying equation (39). Once equation (39) has been 
solved, all flow conditions at station (1) are known and the overall 
continuity equation (34) can be solved. If the calculated flow rate 
is different from the assumed one, V,,) will be changed and the 
iteration shown in Fig. 4 is performed. When overall continuity 
equation is satisfied, the assumed radial position of streamlines 

is checked, and necessary iterations are performed. If streamlines 
continuity is met also, the flow at the rotor inlet is known completely. 
A similar iterative method is then followed for station 2 after the 
rotor as indicated in Fig. 4, If the flow after the rotor is known 
also, the radial shift of the streamlines can be determined, and the 
whole calculating procedure is repeated by taking account of the 
streamlines curvatures and the slopes in solving the equation of 
motion. 

Eventually the resulting average pressure ratio across the turbine 
is compared with the desired one, and if the difference is larger than 
a smallprefixed error, the complete calculating scheme is repeated 
again with a new value of inlet Mach number. After the specified 
pressure ratio has been obtained, all output information is calculated 
and printed. 

FLOW ANGLES AND LOSS COEFFICIENTS CORRELATIONS 

Evidently, even the most accurate method of performance prediction 
will give results in disagreement with actuality, if the flow angles 
and/or the losses in the blade rows are not predicted with the necessary 
accuracy. 

Several correlations, either empirical or theoretical, have been 
proposed in the technical literature. A comparative study of five 
correlations methods available in the present program can be found in 
Reference lk]. ‘The purpose of the present chapter is not to discuss 
the applicability of these different correlations, but only to indicate 
the problems which arose, in applying the different correlation methods. 
In Fig. 5 the nomenclature and sign convection of the velocity triangles 
is indicated. Unless otherwise stated, the variables are positive if 
directed as in Fig. 5. In Fig. 6 the blading geometry characteristics 


used in the different correlation methods are indicated. The names 


eu 


given in Fig. 6 are the FORTRAN names used in the program for stator 
quantities at the mean radius. 

The program is organized in such a way that other than considered 
flow angles and loss coefficients correlations can be introduced by 
adding a few cards and pertinent subroutines. Since the solution method 
is iterative, added subroutines can use blading and flow characteristics 
that are not used in the present program. 

a) Flow angles correlations 

Two methods, one given by Ainley and Mathieson [6] and the other 
given by Traupel [10] are available in subroutines. In both these 
methods,exit flow angles ( 1, Bo )are considered to be independent 
on inlet angles ( Xs» By). While in Traupel's method the angles are 
a function of blading geometry only, in A&M method also the exit Mach 
number is supposed to affect the angles. Although both these methods 
are given for a "reference diameter approach", they will be applied 
in three-dimensional calculations with the assumption that at any 
radius the angle is related only to the corresponding characteristics. 
Since during the calculations both Mach numbers and streamlines 
radii will be changed, flow angles corrections are automatically 
performed during the solution of the equation of motion. 

b) Loss coefficients correlations 

The assumption is made that at any particular radius the loss 
coefficient is a function of the characteristics at that radius only. 
It is realized that this assumption is not realistic, since secondary 
and leakage losses are certainly not equally distributed along the 
radius. On the other hand, secondary and tip leakage flows are not 
understood well enough at the present time to formulate any other 
reasonable assumption. 

Five loss correlation methods given in Ref. [6-10] are available 
in subroutines. They calculate loss coefficients from given blade 
geometry and flow angles and velocities. 

The first considered method is the Ainley and Mathieson |6] 
method. Recently an improvement to the original method was given 
by Dunham and Came 17]. Both methods are available in subroutines. 


These methods express losses by total pressure loss coefficient Y, 
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defined by: 


IZ > IP 
gene Ore se (40) 
P, - Pp 


and assume Y to be independent on pressure ratio Pro/P - The 
relationship between the kinetic energy loss coefficient © and Y is 


given by: 


1+Y dy 


[ 
ace il 
> eee P/Pto = (41) 


The assumption that Y is constant implies that © decreases with increasing 
Mach number. Since this behavior is not experimentally proved, two 
alternatives to equation (41) are available in the program. In both 

these alternatives C is supposed to be constant with pressure ratio, 

and assumes respectively the value corresponding to very low Mach 


numbers (eq. 42) or to an exit Mach number of 0.8 (eq. 43). 


6 = Co = lim 6 = 25 (42) 
Dy Enon L 
vee 
Seah gal 
, See 
( 0 rae (43) 
Ci/pe) Y =1 


where 8* is the pressure ratio corresponding to an isentropic Mach 


number of 0.8, given by: 


a 


px= [1+4*x .87] ae (uy) 
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The same alternatives are given also for the simple method of 
Lenherr and Carter Lo]. 

For the method of Balje' elles the boundary layer wake form 
factor.H must be.introduced. 

With the same nomenclature and assumptions as in eq. (29), H can 


be expressed by: 


Pat je a 
al 
ai (1-X,) il in (222) 
Be FSF (45 ) 
co i-1 co ate 
- (1-X,) 23 ( Xe ) li ye cil Xe ) 
line Tame) 1 L+m(2i-1) 


In this program the assumption m = 0.15 is made. If another exponent 
has to be introduced, card #17 must be changed in the MAIN routine. 
In Traupel's method, the loss coefficients depend on the 


coefficient: 


x = x (Re, relative surface roughness, turbulence factor). (46) 


In the present program, the standard value of a is assumed. If 
the considered turbine has peculiar values of surface roughness or 
Reynolds numbers, it will be necessary to change card #15 in the MAIN 
routine. 

For conformity with the other methods, tip clearance losses are 
expressed by the kinetic energy loss coefficient Core instead of by 
the overall efficiency decrease Allo; as in [10]. The relationship 


between Alloy, and €cz, is given by: 


(47) 


where Ah., is the isentropic enthalpy drop through the turbine and 


Cy is the rotor kinetic energy loss coefficient for zero tip clearance. 
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COMPUTER PROGRAM 


The present computer program is composed of a MAIN routine 
and 28 subroutires, eight of which are concerned with flow calcula- 
tions, two with flow angle correlations, 15 with blade loss 
correlations, and three with numerical interpolations. A brief 
description of the purpose of each routine is given below. 

a) MAIN routine 

1. The primary function of the MAIN routine is to control 
the overall logic of the computer program which is indicated in 
Fig. 4. Furthermore the imput data are read, numerical constants 
are introduced, interpolations concerning stall and blade angles and 
flow areas are performed and the flow angles are computed. Eventually, 
also the output data are printed. 

b) Flow calculation subroutines 

2. In subroutine CHAN the flow rate passing through the section 
O is calculated from the conditions ahead of the stator. 

3. In subroutine STATJR the equation of motion (18a) is solved, 
and the flow conditions after the stator are calculated. If necessary, 
new values of the flow angles are also calculated. 

4, In subroutine ROTORL the absolute flow properties calculated 
in STATOR are converted into relative properties at rotor inlet. 

5. In subroutine RATOR2 the equation of motion (18) is solved, 
and the flow conditions after the stator are calculated. If necessary, 
also new values of flow angles are calculated. 

6. In subroutine FLOWR the flow rate passing through the blade 
sections, either at station 1 or 2, is computed by equation (34), and 
compared with the value calculated in CHAN. If the difference is larger 
than a small prefixed error, a new value of the axial velocity at the 
mean stremaline is computed. The fractional values of the flow rate 
passing between the hub and each streamline are calculated also. 

7. In subroutine SLINE the fractional values of the flow rate 
computed in FLOWR are compared with prefixed ones. If the differences 
are larger than a small specified error, new values of the streamline 


radii, either at station 1 or 2, are computed. 


a 


8, In subroutine ALAS1 the stator loss coefficients are 
calculated from the chosen loss correlation, as a function of 
blading geometry and flow characteristics. Coefficients C* of (28) 
are computed for the stator blades. 

9. In subroutine ALOS2 the rotor loss coefficients are calcu- 
lated from the chosen loss correlation as function of blading 
geometry and flow characteristics. Coefficients ¢* of equation (28) 
are computed for the rotor blades. 

c) Flow Angles Subroutines 

10. In subroutine ANGAIN the flow angles are computed from 
blading geometry and exit Mach numbers with the method of Ainley 
and Mathieson [6]. 

11. In subroutine ANGTRA the flow angles are computed from 
blading geometry with Traupel's method [10]. 


d) Loss Correlation Subroutines 


The six subroutines given below are concerned with the methods 16] and 
har 

ll.» In subroutine COEFFI the imput cards concerning the profile 
losses of nozzle blades given in [6] are read, and the necessary 
interpolations are performed. 

12. In subroutine ChEFF2 a similar process is followed for the 
profile losses of impulse blades. 

13. In subroutine CALYA# the profile losses of nozzle blades are 
calculated with the coefficients computed in COEFFI. 

14. In subroutine CALYA2 the profile losses of impulse blades 
are calculated with the coefficients computed in COEFF? . 

15. In subroutine FIG6 the necessary interpolations for calcu- 
lating the secondary losses are performed. 

16. In subroutine AINLOS the overall loss coefficients with 
either method l6! or Leal are calculated from blading geometry and 
flow angles. Subroutine AINLOS must be called after subroutines 
COEFFI and C(EFF2, and calls subroutines CALYA®, CALYA2 and FIG6. 
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The eight subroutines given below are concerned with method (11) 
of the loss correlation. 

17. In subroutine TRAUP1 the imput cards concerning the 
profile losses given by Traupel's method are read, and the necessary 
interpolations are performed. 

18. In subroutine XP) the basic profile loss coefficients are 
computed with the coefficients calculated in TRAUP1. 

19. In subroutine CSIM Mach number effects on the loss 
coefficients are computed. 

20... In. subroutine CID the effects due to the trailing edge 
thichness and nontwisted blades on the loss coefficients are computed. 

el. In subroutine CSIW the end wall loss coefficients are 
computed. 

22. In subroutine CSIR the secondary loss coefficients are 
computed. 

23. In subroutine ALEAK the tip leakage efficiency losses are 
computed. 

24, In subroutine TRAUP2 the overall loss coefficients are 
computed with Traupel's method Lio] from both blade geometry and flow 
characteristics. Subroutine TRAUP2 has to be called after TRAUP1, 
and calls subroutines XPO, CSIM, CID, CSIW, CSIR and ALEAK. 

25. In subroutine ZA1LOS the overall loss coefficients are 
computed with the method of Lenherr and Carter Lg] from the flow 
angles and the relative velocities. 

26. In subroutine BALJE the overall loss coefficients are 
computed with Balje's method [8] from blade geometry and flow 


conditions. 


e) Numerical Interpolation Subroutines 
27. In subroutine PARAB the coefficients of the parabolic 


interpolations used to approximate curves are calculated. 
28. In subroutine CHBFT the coefficients of the polinomial of 
arbitrary degree, which is the best fit of a given array of points 


in the Chebyschev sense, are calculated. 


if 


29. In subroutine YC the value of the polynomial determined 
in CHBFT, corresponding to a given value of the independent 
variable, is calculated. 
Description of Imput Data 

The imput data of the present computer program consist of 63 + 
NSETS cards, NSETS being the number of conditions for which 
performance values are to be obtained. They can be divided in three 
sections. In the first, the constant values used in the methods 
of Ainley and Mathieson and of Traupel are introduced (these 36 
cards are read respectively in subroutines COEFFI, CCEFF2 amd TRAUP1). 
In the second, the specifications about the detail of output desired, 
the tolerances of the numerical solutions, the systems of correlatim 
of the flow angles and the loss coefficients, the coefficients used in 
the equation of continuity, the gas properties am the turbine geometry 
are stated. In the third, the conditions for which performance values 
are to be found are introduced. 

The 36 input cards for the first section are given in Table I. 
They will remain unchanged for all turbines. Detailed description 
of the 27 cards composing the second section of the imput are given 
in Table II, where also a description of the third imput section can be 
found in cards #28 and #29. An example of the imput cards described 
in Table II is given in Table III. 
Description of the normal output 

In Table IV the output corresponding to the imput of Tables I 
and III is given. It can be seen that it consists of 1+ 3 NSETS 
sheets. In the first sheet the turbine geometrical dimensions and 
the choice among the different correlation methods are described. 
Then three sheets for every set of operating conditions are printed. 
The first of these sheets describes the flow conditions at station 1 


after the stator row; the second one the flow conditions after the 
rotor row; and the third one gives the main overall turbine characteristics 


and the mass averaged performance values. 
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Fig.4 Flow diagram 
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CORRELATION POSSIBLE 


SPECIFICATIONS 
INDEX | VALUES 
-= be fe Ne a ae 
oe | Asm [6] or D&C [7] loss correlation | 
| method is adopted, depending on IAI | 
| value. 
| 
2 Balje [8] loss correlation method is © 
| adopted. | 
| | 
: 3 L&C [9] loss correlation method is : 
! adopted. 
| 
h Traupel [10] loss correlation method 
| is adopted. 
a0 _ If IcOR=1, A&M [6] loss correlation 
| method is adopted. 
i If ICOR=1, D&C L7] loss correlation 
method is adopted. 
il | Asm L6] flow angles correlation is 
adopted. 
IAN 
| 2 Traupel Lio] flow angles correlation 
| is adopted. 
—————— in ie a O 
| =o Rotor is not shrouded. 
| | 
aor | 
ak Rotor is shrouded. 
H 
TABLE V (1) 


yy 


CORRELATION 
INDEX 


POSSIBLE 


SPECIFICATIONS 


Loss coefficients are function of 


flow angles 
Loss coefficients are function of 


blade inlet angles, for negative 


incidence 


Loss coefficient C = C, (see eq. (42)) 


Loss coefficient ¢ = C (Y, pressure 
ratio) (see eq. 41)) 


Loss coefficient ¢ = ew-0.8 (see eq. (43)) 


coefficient in continuity equation 


-2 Cror 


coefficient in continuity equation 


CPROF 


coefficient in continuity equation 


Cron 


TABLE V (2) 
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